Oscillatory athermal quasi-static deformation of a model glass 
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We report computer simulations of oscillatory athermal quasi-static shear deformation of dense 
amorphous samples of a three dimensional model glass former. A dynamical transition is observed 
as the amplitude of the deformation is varied: for large values of the amplitude the system exhibits 
diffusive behavior and loss of memory of the initial conditions, whereas localization is observed 
for small amplitudes. Our results suggest that the same kind of transition found in driven colloidal 
systems is present in the case of amorphous solids (e.g. metallic glasses). The onset of the transition 
is shown to be related to the onset of energy dissipation. Shear banding is observed for large system 
sizes, without, however, affecting qualitative aspects of the transition. 



Understanding the behavior of materials under me- 
chanical deformation is of primary importance for many 
contexts. While the deformation behavior of crystals 
is theoretically well understood, no universally accepted 
framework exists to rationalize the behavior of mechan- 
ically driven amorphous systems, although significant 
progress has been made in recent years in developing a de- 
tailed understanding of how an amorphous solid responds 
to external stresses [TH3]- Considerable recent activity 
has been spurred by an interest in the mechanical behav- 
ior of metallic glasses, soft glassy materials, foams and 
granular packings, and has involved theoretical, compu- 
tational and experimental investigations [21 IMS] ■ Partic- 
ular interest is understandably focused on the manner in 
which the response of an amorphous solid changes from 
nearly elastic response at small applied stress to a state 
of flow for large applied stress. 

Many computational investigations have employed the 
approach of studying the zero temperature behavior of 
amorphous solids under quasi static conditions (using an 
Athermal Quasi Static or AQS procedure [S]). In this 
procedure, systems are kept in local energy minimum 
configurations, or inherent structures [101 E] while vary- 
ing the strain. In previous work on model systems of 
binary Lennard- Jones particles, it has been shown that 
upon monotonically increasing the applied shear strain, 
the inherent structures evolve towards energies corre- 



sponding to the limit of high temperatures [T^] . This and 
related phenomena are referred to as rejuvenation, in con- 
trast to the well studied process of ageing whereby (typi- 
cally) a glassy material descends to deeper energy config- 
urations as a function of the waiting time over which it 
relaxes at a given temperature. In contrast, when a cy- 
cle of strain is applied up to a maximum value which is 
then reversed |13j . both ageing and rejuvenation are ob- 
served, with small amplitude strains found to reduce the 
energy of samples ("overage" them), while larger ampli- 
tude strains tend more often to increase the energy (thus 
"rejuvenating" the samples). Initial conditions of the 
samples in such cases matter: samples with lower initial 
potential energy are rejuvenated more easily than those 
with higher energy I3J. 

In a very different context, it has been observed that 
the nature of rearrangements of colloidal particles in a 
suspension, at moderate packing fractions, subjected to 
cyclic shear deformation, depend on the amplitude of 
strain [T3]. For small amplitudes the particles become 
quiescent after a short transient, whereas for larger am- 
plitudes a finite fraction of the colloidal particles move 
irreversibly [HI [IS] . It has been argued that such a non- 
equilibrium transition from a quiescent to an active state 
belongs to the universality class of conserved directed 
percolation (C-DP) [TF. 

In this letter we report results of computer simu- 
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lations for a three dimensional binary Lennard- Jones 
model, wherein we subject samples to a large number 
of shear strain cycles (rather than to a single cycle as 
reported in |13p. Our systems reach, for large enough 
amplitudes, steady states which are independent of 
their initial conditions. Interestingly, initial systems 
populating different regions of the potential energy 
landscape converge to the same steady states with an 
energy that depends only on the maximum strain. In 
such cases the systems exhibit diffusive behavior, with a 
diffusion coefficient that is characteristic of the steady 
state, depending only on the strain amplitude. For low 
amplitudes of the strain, samples with different initial 
conditions fail to converge to a common steady state, 
and the dynamics is not diffusive, but samples reach 
instead non-diffusive, quiescent states that strongly 
retain memory of the initial conditions. Our results 
thus indicate that sheared amorphous solids undergo 
the same type of quiescent to active transition as driven 
colloidal suspensions at moderate density [15 . We find 
that particle diffusivity is a suitable order parameter 
for such transition, and suggest that the same diffusion 
behavior observed in driven colloidal suspensions could 
occur in completely different cases, namely in fatigue 
experiments on metallic glasses [5l[T7l[T8]. In addition, 
we find that the critical amplitude 7c at which the 
transition is observed depends on the size of the samples 
studied, and is directly related to the dramatic change in 
the amount of energy dissipated per deformation cycle 
(the area of the hysteresis curve), rather than to the 
value of the yield strain of undeformed samples. Finally, 
we show that large systems subjected to deformation 
above the critical value present indication of shear 
banding. Remarkably, while the band can occupy as 
much as half of the simulation box, the above dynamical 
transition picture remains unaffected. 

We consider a binary Lennard- Jones mixture with 
parameters as in fl9 . Samples with N = 500, 4000, 
32000 particles are equilibrated at two different tem- 
peratures (T — 1.0,0.466) at reduced density 1.2. The 
higher temperature corresponds to the onset of slow 
dynamics whereas the lower temperature, close to the 
lowest temperatures studied for this system, corresponds 
to a supercooled state [20]. Configurations sampled 
from the equilibrium trajectory are subjected to energy 
minimization, using the conjugate-gradient algorithm, to 
obtain sets of inherent structures typical of the liquid at 
temperature T. The inherent structure energies sampled 
by the liquid are lower for the lower temperatures |13j . 
All the simulations were performed using LAMMPS [21] . 

Each of such inherent structures is then subjected to 
two operations: first, particle positions are transformed 
affinely (in this case, we transform the x-coordinates ac- 
cording to = + d'-fTy while keeping Vy and the 



same) so that the shear strain 7 of the samples is incre- 
mented by a small amount dj, and Lees-Edwards bound- 
ary conditions [22] are updated accordingly. This has the 
consequence of changing the original energy landscape 
in which the inherent structure configuration is not any 
longer a minimum. Subsequently, energy minimisation is 
performed using the conjugate gradient algorithm, to ob- 
tain the inherent structure in the new, deformed, energy 
landscape. These two operations represent an Ather- 
mal Quasi-Static (AQS) step [^. If AQS steps are it- 
erated, the sample can be sheared to arbitrary values 
of strain. The dynamics under such a protocol can be 
assumed to mimic real shear deformation experiments 
when the temperature of the system is low enough to 
be neglected (hence the name athermal) and shear rates 
are low enough to let the system relax fully before fur- 
ther deformation takes place (so that the deformation 
can be considered quasi-static) [9]. In what follows, the 
shear strain 7 is changed in steps of d'y {— 10^^, 2 • 10~*, 
2 • 10-"*, respectively for N = 500, 4000, 32000) from to 
a maximum value ^max and then to —jmax and back to 
0. This periodic operation is the elementary cycle of our 
dynamics. For practical purposes we define the accumu- 
lated strain 7acc as the sum of the elementary absolute 
strains applied at all steps: jacc — J2i MTjI- The accu- 
mulated strain plays the role of time in our analysis. At 
the end of each cycle, particle positions in configurations 
with 7 = (zero strain configurations) arc stored, so that 
the mean squared displacement with respect to the initial 
configuration (7acc — 0) can be measured. Energy and 
all the components of the stress tensor are also recorded 



at every AQS step (but, except in 3b we report only the 
zero strain values at the end of each cycle). Up to 40 
samples are subjected to cyclic shear deformation using 
the protocol above, so that the quantities above can be 
averaged. 

The average potential energy of the minima visited 
is presented in Fig. [Taj as a function of the number of 
cycles. The behavior of the potential energy depends on 
the value of jmax- For low values of jmax, after a short 
transient, the energy curves reach a value which depends 
on the initial T. For high enough values of jmax, the 
energy converges to the same 7mas-dependent value 
for both the initial values of T, and remains steady. 
We expect this result to hold for other temperatures 
as well. Whether overaging or rejuvenation will occur, 
therefore, depends on whether the initial energy of 
the undeformed inherent structure is lower or higher 
than the steady state energy. At intermediate values of 
Jmax, the convergence of the energy to the asymptotic 
values requires more cycles than at the extreme jmax 
values. To characterize this "slow relaxation" , we fit the 
transient behavior of the energy data in Fig. [Taj with 
a stretched exponential, from which a characteristic 
relaxation strain "face is extracted. As can be seen from 
Fig. |lb[ the values of jacc become high for intermediate 
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Imax and show behavior consistent with a divergence at 
a critical strain, similarly to what is observed in driven 
colloidal suspensions [15]. 
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FIG. 2. Mean squared displacement for configurations, for 
different '^rnax, averaged over different runs with A'^ = 4000 
particles at T = 0.466 (filled symbols) and 1.0 (open sym- 
bols) (Symbols as in Fig. lb I. Note the transition between an 
arrested and diffusive regime as "fmax is increased. 
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FIG. 1. a) Potential energy per particle E for zero strain con- 
figurations, for different values of 'ymax ((O)-0-07, (0)10. 08, 
(O):0.09 (V):0.1 (A):0.12 (*):0.14), averaged over diflferent 
runs with samples of iV = 4000 particles at T = 0.466 (filled 
symbols) and 1.0 (open symbols). The fits are obtained us- 
ing a stretched exponential model. For the highest values of 
7maa:, E fluctuates arouud a plateau after an initial transient, 
whose value increases with "fmax and is not dependent on the 
starting energy (curves with filled and open symbols of the 
same color/symbol merge), whereas for lower "fmax memory 
of initial conditions is retained. For intermediate values our 
simulations are not long enough to determine which of the two 
behaviors eventually holds, b) Values of characteristic strain 
7acc (akin to relaxation time) obtained from the stretched ex- 
ponential fits of the energy relaxation in [Ta| The relaxation 
strain increases as one approaches the transition between the 
arrested and diffusive regimes. The lines through the data 
are guides to the eye. Vertical lines indicate critical strain 
amplitudes 7c described in Fig. |3a[ 

In order to understand the nature of particle motions 
that lead to the energy relaxation (or lack thereof) that 
we observe, we calculate the mean squared displacements 
(MSD) from the initial configurations, averaged over the 
samples. The MSD curves shown in Fig. [2] indicate that 
for small values of ^max the MSD curves saturate after an 
initial transient, whereas for large ^max the particle dis- 



placements are diffusive, i.e., the MSD depends linearly 
on the accumulated strain beyond some value of 7acc- In 
the case where a steady state is reached, the diffusivity 
can be extracted by a linear fit of the MSD calculated 
from an initial point in the steady state. In terms of 
the relative accumulated strain to such a point, 7*f,c, the 
MSD is given by 



MSD = D-fl 



(1) 



At intermediate values of ^max^ the system does not 
reach either an immobile or a diffusive steady state within 
the number of cycles we perform. The energy in Fig. la 
and the MSD data in Fig.|2j taken together, suggest that, 
in this regime, the system undergoes a transition from 
an arrested (or quiescent) to diffusive behavior across a 
critical ^max, near which the relaxation of the system 
becomes very sluggish. 



In Fig. 3a the values of D obtained for all the system 



sizes are plotted against jmax- For low 'fmax, D is zero 
and starts to increase rapidly around a critical value 7^ 
which depends on the system size (Fig. [3a]). The values of 
D for ^max > 7c in Fig. [3a] are reasonably well described 



by a law of the form (as in \T5\ ) 



Ic) 



(2) 



The relevant fit values are mentioned in the caption of 
Fig. [3a[ Although better analysis is needed to be sure of 
the values of the exponent /3 and critical strain 7c, the 
data shown are clearly consistent with a transition from a 
regime with zero diffusivity to one with finite diffusivity. 
In the large strain regime, particles move diffusively and 
explore configuration space whereas, in the low strain 
regime, they are localized in configuration space. To test 
the robustness of the observation of a finite 7c, we plot 
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7c for different system sizes in the inset of Fig. 3a against 
1/N . This plot indicates that 7c will be finite for N oo. 
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FIG. 3. a) Values of diffusivity D obtained by fitting the 
MSD, for the studied system sizes and temperatures. Super- 
imposed on the data are power law fits (Eq. [2| from which the 
critical strain (7c) values obtained (for A'' = 500, 4000, 3200 
respectively, for (T = 0.466, T = 1.0)) are: (0.115,0.117), 
(0.083, 0.083) and (0.070, 0.073). The corresponding exponent 
values P are: (0.38,0.36), (0.75,0.76), (0.61,0.54). The inset 
shows 7c vs. inverse system size 1/N. 7c extrapolates to a 
non-zero value as 1/N — >■ 0. b) Stress-strain curves (shown in 
the inset for N = 4000, T = 0.466) show hysteresis with finite 
enclosed area that increases with ^max at high strain ampli- 
tude "fmax ■ Also shown for reference is the stress-strain curve 
obtained under uniform shear. The areas enclosed by the hys- 
teresis curves are shown in the main panel as a function of 
"fmax for all the different system sizes and temperatures stud- 
ied. The areas, close to zero at low fmax, become finite above 
7c (indicated by vertical lines for each system size). 



Interestingly, the transition to the diffusive regime also 
coincides with the onset of energy dissipation in the sys- 
tem, as can be deduced by considering the stress-strain 
curves in the steady state, which are shown in the in- 
set of Fig. 3b For small jmax, these curves are nearly 



linear, and do not enclose any significant area, indicat- 
ing nearly elastic behavior. For large strain however, the 
curves begin to clearly show hysteresis. This effect is 
common in strain experiments on several materials rang- 
ing from metals to bio-materials and it is an indication 



of a viscoelastic behavior. As shown in Fig. 3b the ^max 
values beyond which dissipation becomes finite agree, to 
a very good extent, with the critical strain 7c for the 
corresponding system sizes. These results therefore show 
that the transition from a localized to a diffusive regime 
corresponds also to a transition from an nearly-elastic 
(almost non-dissipative) to a dissipative regime. It is in- 
teresting to ask if the critical strain that we observe is 
related to the yield strain in steady shearing conditions. 
Although such a possibility is supported by results for 
the smallest system we study, the yield strain (which we 
find not to be very sensitive to system size) is higher than 
the critical strain for the larger ones {N — 4000, 32000) 
(see inset of Fig. 3b I. 

Although our results for different system sizes show 
qualitatively consistent behavior, at the largest size stud- 
ied, N = 32000, a new phenomenon is observed. The 
system under strain displays clear indications of shear 
banding. Fig. |4] shows a snapshot in which particles that 
have moved by more than 0.6a in one cycle are high- 
lighted. Such particles are clearly correlated in position, 
forming a shear band. Shear bands are present in all 
the samples above the critical strain. However, despite 
the emergence of shear banding, the dynamical transition 
scenario discussed so far is not affected. 




FIG. 4. Snapshot of a A'^ = 32000 configuration obtained in 
the steady state for 7max = 0.08. Particles in red are those 
that have undergone a scalar displacement > 0.6a when the 
sample is subject to a full shear deformation cycle on the xy 
plane. The presence of a band is evident. 

In conclusion, we have shown that athermal oscillatory 
shear deformation drives, for large enough strains, dense 
samples of a model glass into steady states in which they 
explore their energy landscape in a manner dictated by 
the deformation amplitude jmax- Previous observations 
of overaging and rejuvenation via a single cycle of de- 
formation in a model glass |13) are thus rationalized as 
initial steps towards a steady state which is independent 
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of the initial state of the samples and depending only on 
Imax- Depending on whether the amplitude of the de- 
formation is below or above a critical strain amplitude 
7c, these driven amorphous solids can be in an arrested 
(or localized) state or in a diffusive state which is char- 
acterised by hysteresis in the steady state stress-strain 
curves. The relaxation to the steady state becomes very 
sluggish for ^max values near the critical value. The tran- 
sition from a localized to a diffusive steady state we ob- 
serve is very reminiscent of that observed in experiments 
on driven systems [151 123] and reported recently by other 
authors on systems similar to ours. The value of 

the amplitude 7c at which this transition occurs depends 
on the system size but appears to remain finite in the 
limit of infinitely large systems. We find that the value 
of 7c has an excellent agreement with the strain value at 
which hysteresis develops in the stress-strain curves, indi- 
cating a transition from a quasi-elastic (non-dissipative) 
to a viscoelastic (dissipative) regime. The critical strain 
values, however, do not coincide with the yield strain, 
even though these values are close. Shear banding oc- 
curs for our largest system sizes {N > 32000) in diffusing 
samples. Interestingly enough, the onset of shear band- 
ing at large N doesn't fundamentally affect the dynami- 
cal transition. 
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